Dilute Bose gases interacting via power-law potentials 
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Neutral atoms interact through a van der Waals potential which asymptotically falls off as r -6 . 
In ultracold gases, this interaction can be described to a good approximation by the atom-atom 
scattering length. However, corrections arise that depend on the characteristic length of the van der 
Waals potential. We parameterize these corrections by analyzing the energies of two- and few-atom 
systems under external harmonic confinement, obtained by numerically and analytically solving the 
Schrodinger equation. We generalize our results to particles interacting through a longer-ranged 
potential which asymptotically falls off as r -4 . 

PACS numbers: 



I. INTRODUCTION 



The interaction strengths of sufficiently dilute and cold 
bosonic atom samples such as Bose-Einstein condensates 
of alkali atoms can be parameterized to a good approx- 
iraation by a single parameter, the s-wave scattering 
length 1]. In these systems, the neutral atoms inter- 
act through short-ranged van der Waals potentials which 
fall off as r~ 6 at large interparticle distances r. More 
recently, progress has been made in cooling and trapping 
systems characterized by interaction potentials that fall 
off more slowly than r -6 at large r. For example, the 
interaction between a neutral atom and an ion is domi- 
nated by a polarization potential that falls off asymptot- 
ically as r -4 0. Atom-ion systems have recently been 
proposed as candidates for quantum computing applica- 
tions and also play a role in recent work which pro- 
poses that macroscopic molecules can be formed by im- 
mersing an ion in a condensed Bose gas 0,(1,0. Another 
example for systems with longer-ranged interactions are 
dipolar gases [7|, |8|. In these systems, the non-neglegible 
magnetic or electric dipole moment leads to an angle- 
dependent r~ 3 potential at large interparticle distances. 
A natural question to ask is how well the properties of 
Bose systems with longer-ranged interactions can be de- 
scribed by the s-wave scattering length. 

This paper considers dilute bosonic systems under ex- 
ternal confinement interacting through spherically sym- 
metric power-law potentials. In particular, we treat in- 
teractions with r~ n tails, where n is 4 or 6. We fo- 
cus on the regime where the characteristic length (3 n of 
the two-body potential is much smaller than the char- 
acteristic length ciho of the trapping potential. In this 
regime, the shape-dependent interaction potential can be 
replaced by a regularized zero-range potential whose in- 
teraction strength is parametrized by the s-wave scat- 
tering length. For the potential with r -6 tail, e.g., it 
has been shown previously that the energy levels of the 
trapped two-body system can be reproduced very accu- 
rately if the energy-dependence of the scattering length 
is accounted for [9|, [lfj]. This paper extends the two- 
body analysis to potentials with r~ 4 tail, whose scatter- 
ing length has — because of the longer-ranged character of 



the potential — a stronger energy-dependence than that of 
potentials with r~ 6 tail. We find that the corrections to 
the energy predicted by the zero-energy scattering length 
go as (Pe/aho) 3 and (/^/a^o) 2 for the interaction poten- 
tials with r~ e and r -4 tails, respectively. 

Using Monte Carlo techniques, we furthermore treat 
dilute bosonic many-body systems. As in the two-body 
case, we consider different interaction potentials and an- 
alyze the resulting eigenenergies. Not unexpectedly, our 
results show that the energy-dependent scattering length 
remains a good quantity also in the many-body system. 
This suggests, e.g., that the description of dilute Bose 
gases within a mean-field Gross-Pitaevskii framework can 
be improved notably by including the energy-dependence 
of the scattering length. First steps in this direction have 
already been taken [y, [ll[ ; our results provide additional 
benchmark results that may aid in further assessing the 
accuracy of these and related frameworks. 

Section |TT] introduces the Hamiltonian and the model 
interaction potentials used in our study. Section III dis- 
cusses the energetics of two particles in a trap interacting 
through both finite-range and zero-range potentials. In 
Sec. IV, we consider the energetics of more than two par- 
ticles in a trap by solving the many-body Schrodinger 
equation using Monte Carlo techniques. Finally, Sec. V 
concludes. 



II. HAMILTONIAN 

The Hamiltonian for a system consisting of N identical 
mass m bosons in the presence of a spherically symmetric 
harmonic trapping potential with angular frequency to is 
given by 
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where denotes the position vector of the ith atom. The 
spherically symmetric two-body interaction potential v 
depends on the relative distance ry, r,j = |rj — Tj\, We 
consider attractive power-law potentials with a hardcore 
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with n = 4, 6 and C„ > 0. The Hamiltonian defined 
in Eq. (|TJ) is characterized by three length scales: the 
hardcore radius r c , the van der Waals length scale (3 n 
{Pn — (mC n /fi 2 ) 1 '''"" 2 '], and the harmonic oscillator 
length a/jo [a^ = y/h/muj]. Throughout this paper, we 
are interested in the regime where r c and (3 n are much 
smaller than a^ Q . 

In three dimensions, the interaction strength of any 
potential that falls off faster than 7 ,_3 at large distances 
can be characterized by the e nerg y-dependent free-space 
s-wave scattering length a(k) 
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where S(k) denotes the s-wave scattering phase shift and 
k the wave vector at the scattering energy of H 2 k 2 /m. 
The zero-energy scattering length a(0) is defined by tak- 
ing the k — > limit of Eq. ([3]). For the V4 and vq po- 
tentials, the zero-energy and energy-dependent scatter- 
ing lengths can be calculated from analytical solutions 
derived using series expansion techniques [l2l . [l3| . 

Figures HJa) and [21a) show the zero-energy scatter- 
ing length a(0) as a function of [3 n for the v n potential 
with r c = 0.007aho for n equals 4 and 6, respectively. 
Although the scattering lengths are calculated for the 
free-space system with no external trapping potential, we 
choose to express all lengths in units of a^o to ease the 
comparison with the trapped system in Sees. Ill and IV. 
When (3 n = 0, the scattering length coincides with the 
hardcore radius r c . As j3 n increases, the attractive tails 
of the v n potentials increase in strengths, which leads to 
a decrease of the scattering lengths. This continues until 
the potential is strong enough to support its first bound 
state, at which point the scattering length changes its 
sign from negative to positive. Figures [TJa) and[2ja) in- 
dicate that this divergence occurs at different values of 
/?„, i.e., at (ii ~ 0.022a/ lo and /3 6 « 0.016a^ o , owing to 
the fact that the v§ potential is shorter-ranged than the 
t>4 potential. This can be understood heuristically by 
considering the ratio of the attractive power of the po- 
tentials, that is, the ratio J ^6(r)d 3 r/ J V4(r)d 3 r. Taking 
the limits of the r-integration as r c and 00, this ratio 
equals /3 4 /(3/3f r;?). Looking at equal values of (3^ and (3q 
and considering that (3 n > r c at the first divergence, this 
ratio is greater than one, in agreement with the obser- 
vation that the vq potential supports an s-wave bound 
state for smaller values of (3 n /aho than the V4 potential. 

Throughout this paper we are interested in describ- 
ing dilute Bose systems which interact primarily through 
binary s-wave collisions. In such systems, the short- 
range details of the interaction potential are not being 
probed, and the regularized zero-range pseudopotential 



reproduces many observables obtained for the true shape- 
dependent interaction potential — in our case, the V4 or vq 
potential — accurately if the strength g is chosen properly. 
In Sees. Ill and IV we take g to be the zero-energy scat- 
tering length a(0) and the energy-dependent scattering 
length a{k) of the shape-dependent interaction potential 

Vn- 

For two particles in a harmonic trap interacting 
through v ps , the Schrodinger equation can be solved 
analytically [l5j ]. The center of mass energy equals 
(n cm + 3/2)hu> (n cm = 0, 1,. . .), and the s-wave eigenen- 
ergies E re i = efno of the Schrodinger equation in the 
relative coordinate are determined by [15j 
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The transcendental equation (JSJ) can be solved straight- 
forwardly for any given g using standard root-finding 
procedures. For N > 2, analytical solutions to 
the Schrodinger equation for trapped atoms interacting 
through v ps are in general not known, and we instead 
resort to numerical techniques (see Sec. IV). 



III. TWO PARTICLES IN A TRAP 

We first consider the Hamiltonian given in Eq. (fTJ) with 
V = v n for N — 2. After separating off the center of 
mass motion, we are left with a Schrodinger equation 
in the relative coordinate. We solve the corresponding 
one-dimensional differential equation numerically using 
£?-splines. Figures [2(b) and[2jb) show the resulting rel- 
ative s-wave eigenenergies, denoted by e V4 and e Vfs , re- 
spectively, as a function of (3 n - As in Figs, [lja) and 
[5Ja), the hardcore radius is fixed at r c = 0.007a^ o . For 
those (3 n values for which a(0) is small [see Figs. Ufa) 
and UK a)], the eigenenergies e Vn coincide approximately 
with the eigenenergies (2n re ;+3/2) of the non-interacting 
system, where n re i — 0, 1,. . .. However, each time a(0) 
diverges, a new molecular state appears and the energy 
of the gas-like state decreases by approximately 2hco. 

Next, we consider two trapped particles interacting 
through v ps . We find that the zero-range pseudopoten- 
tial with energy-dependent scattering length reproduces 
the eigenenergies e Vn for the shape-dependent potential 
v n with high accuracy for all interaction strengths con- 
sidered in Figs. [1] and O To obtain the eigenenergies for 
v ps with g — a(fc), which we denote by e a ( e ), we calculate 
a(k) for different (3 n and solve Eq. §5§ self-consistently 
S LkLHl i i-e., we require that efku on the right-hand side of 
Eq. (O agrees with the energy h 2 k 2 /m at which the two 
particles collide. Since e Vn and e a ( £ ) coincide to many 
digits, Eq. © can be used to describe the physics of 
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FIG. 1: (Color online) s-wave properties of two particles in- 
teracting through the potential V4 with r c = 0.007(^0 as a 
function of j3i/aho (note that the z-axis is the same for all 
three panels): (a) Free-space zero-energy scattering length 
a(0). (b) Relative energy e„ 4 for two trapped atoms, (c) En- 
ergy difference Ae, Ae = e„ 4 — £ a (o)> f° r two trapped atoms. 
In (b) and (c), the line styles are keyed to each other for ease 
of comparison. 



FIG. 2: (Color online) s-wave properties of two particles in- 
teracting through the potential v§ with r c = 0.007ah o as a 
function of Ps/a^o (note that the x-axis is the same for all 
three panels): (a) Free-space zero-energy scattering length 
a(0). (b) Relative energy e V6 for two trapped atoms, (c) En- 
ergy difference Ae, Ae = e V6 — £ a (o)i f° r two trapped atoms. 
The inset of (c) plots Ae on an enlarged scale to show the 
maximum of Ae. In (b) and (c), the lines styles are keyed to 
each other for ease of comparison. 



two trapped particles provided j3 n <C a; lo and provided 
the energy-dependence of the scattering length is known. 
For short-range potentials, this was already shown in 
Refs. 

For some systems under experimental study, 
only a(0) is known [a(k) is unknown]. It is thus useful 
to quantify the deviations Ae between the eigenenergies 
e Vn and the eigenenergies £ a (o) obtained from Eq. (JSJ) with 
g = a(0). 

Figures [He) and \2lc) show the energy difference Ae, 
Ae = e Vn — e a ( ) , for the three energetically lowest-lying 
gas- like states. The line styles in Figs.QTc) and[2jc) cor- 
respond to those used in Figs.QJb) and^b). The energy 
difference Ae is larger for the energetically higher-lying 
gas-like states since the energy-dependence of a(k) for a 
given P n increases with increasing e. The maximum of 
Ae increases with increasing (3 n /aho [see Fig. QJc) and 
the inset of Fig. \2I.c)] and is of the same order of magni- 
tude for the «4 and vq potentials for comparable values of 



Pn/aiio- Furthermore, for those f3 n values for which the 
scattering length a(0) is comparatively small, the energy 
difference Ae also increases with increasing (3 n /aho- The 
magnitude of these "background" energy differences is 
much larger for the potential than for the vq potential 
[note the difference in the y-scales of Figs. QTc) and^c)]. 
For example, for f3 n ~ 0.03a/ lo , the background energy is 
about 10~ 3 hoj and 10~ 5 ?mj for the lowest-lying gas-like 
states of the 1*4 and vq potentials, respectively. We now 
show that the background energy difference Ae is pro- 
portional to {fti/aho) 2 and {(3s/aho) 3 for the V4 and v§ 
potentials, respectively, thus explaining the much smaller 
energy difference for the v§ potential than for the V4 po- 
tential. 

To arrive at these estimates, we use that £ a (e) agrees 
to many digits with e Vn , which implies Ae = £ a (e) ~ £ a(o)- 
Since £ a (e) is determined from Eq. (JSJ) with g — a(e), we 
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FIG. 3: (Color online) Energy difference As, As — e v „ — e a (o), 
for the three energetically lowest-lying gas-like states of two 
trapped atoms for a(0) = as a function of f3 n /a>ho- The filled 
(open) triangles, circles and diamonds show the numerically 
determined energy differences for the levels near 1.5, 3.5 and 
5.5/kJ, respectively, for the v§ (vi) potential. The solid and 
dashed lines show the corresponding analytically determined 
estimates for Ae, Eqs. (fT3)l and ([12}. 



can obtain a simple expression for £ a (e) by expanding the 
left-hand side of Eq. ([5]) about a(0) and the right-hand 
side about £ a m). The expansions of a(k) for the and 



the vq potential are given by [If 



and m 



a{k) = a (0) + -/3 2 fc + 



a(k) = a(0) ( 1 + -r e a(0)fc 2 + 
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respectively. In Eq Q, r e denotes the effective range of 
the vq potential [ll|, [l8| , 
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where the constant x e = [r(l/4)] 2 /(2?r) w 2.09. Denot- 
ing the right-hand side of Eq. |5]) by /(e Q (o)) for = a(0), 
we find 
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for the «6 potential. The different powers of e a (o) 
in Eqs. © and (flT))) follow directly from the linear 
and quadratic /c-dependence of the correction terms in 
Eqs. ([6]) and ((7]), respectively. 

If a(0) <C /?6i the square bracket in the expression for 
the effective range r e is approximately equal to 2 and 



Eq. (|T0[) reduces to 
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Furthermore, for small <z(0), e a (o) is approximately given 
by (3/2 + 2n re i) and /'(e (o)) is of order 1 (taking values 
of approximately 1.25, 0.84, and 0.67 for n re i = 0, 1, and 
2). The second term in the denominator of Eqs. © and 
(fTTj) can thus be dropped provided f3 n is much smaller 
than aho- This yields 



_ iiyj2n rel + 3/2 / /3 4 
£W 3/'(2n rei + 3/2)Va fto 



for the i>4 potential and 
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(13) 



for the «6 potential. Equations (fl"2"|) and (|13p can also 
be derived by applying first order perturbation theory 
to the trapped two-body system interacting through a 
zero-range potential (see Sec. IV). 

Solid and dotted lines in Fig. [3] show the energy differ- 
ences Ae predicted by Eqs. (fT3"j) and (fT2"|) for n re i = 0, 1 
and 2 (from bottom to top) as a function of [3 n / 1 a^o for 
the Vq and u 4 potential, respectively. For comparison, 
filled and open symbols in Fig. [3] show the correspond- 
ing numerically determined energy differences Ae for the 
three energetically lowest-lying gas-like states. To cal- 
culate these energy differences, we fix r c and j3 n so that 
a(0) = 0, and vary the harmonic oscillator length. We 
find that, as long as r c and j3 n -C dhoi the results shown in 
Fig. [3] are independent of the number of bound states sup- 
ported by the shape-dependent power-law potential v n . 
Figure [3] illustrates that the estimates given in Eqs. fT2")) 
and (fT3)) are quite accurate for small a(0). Thus, our 
derivation shows that the different powers of the char- 
acteristic length scale (3 n , which explain the larger back- 
ground values of Ae for the V4 potential compared to the 
vq potential, can be traced back to the different energy 
dependence of a(k) for the W4 and vq potential. 



IV. N PARTICLES IN A TRAP 

To solve the time-independent Schrodinger equation 
for more than N ~ 2 trapped particles, we resort to 
the variational Monte Carlo (VMC) and diffusion Monte 
Carlo (DMC) techniques 

In the VMC method, the variational many-body wave 
function i\)y is written in terms of a set of variational 
parameters p, which are optimized so as to minimize the 
variational energy Ey, Ey = (ipv\H\ipy) / (ipy\ipy) . The 
energy expectation value Ey is calculated for a given p 
using Metropolis sampling. Motivated by the structure 
of the Hamiltonian H, Eq. (Q]), we write tpy as a product 
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of one-body terms ip and two-body Jastrow terms F [2C 
ii, 



A' 



Vv(ri,...,rjy) = H^n) J]>(r y ), 



(14) 
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where <p(r) = exp(— pir P2 ). The functional form of the 
two-body Jastrow factor F is motivated by the functional 
form of the interaction potential v and by the fact that 
we are interested in describing the energetically lowest- 
lying gas-like state of the many-body Hamiltonian. We 
use 



F(r) = 



fl 



for v 
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for 

F(r) = i sin(fcr + d) for 
ei + e 2 exp(-p 6 r) for 



r < c 

c < r < p 5 

r > p 5 
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for v = v ps . The parameters b, c, d and k are chosen so 
that ipv obeys the boundary conditions implied by the 
many-body Hamiltonian H and so that ipv has the de- 
sired symmetry (see below), while the parameters e\ and 
e 2 are determined by requiring that F and its derivative 
are continuous at r = p§. For each interaction poten- 
tial v, we optimize the variational parameters p, i.e., p\ 
through pi for v = v n , and pi, P2, ps and pe for v = v ps . 

To go beyond the variational calculations, we apply 
two different variants of the DMC algorithm both of 
which use the optimized wave function ipy as a guiding 
function [li|. When the many-body Hamiltonian does 
not support any states with negative energy, the lowest 
lying gas-like state coincides with the true ground state 
of the system. In this case, ipv is nodeless and the DMC 
algorithm with importance sampling results in the ex- 
act many-body energy. When the many-body Hamilto- 
nian supports negative energy states, i.e., molecular-like 
bound states, the energetically lowest lying gas-like state 
possesses nodes, which are imposed in the DMC method 
with importance sampling through the variational wave 
function ipv- This DMC variant, referred to as fixed-node 
DMC method [H, [23[ , determines the lowest energy of a 
state that has the same symmetry as ipv- Importantly, 
the FN-DMC energy provides an upper bound to the ex- 
act eigenener gy o f the excited gas-like state of the many- 
body system [231 ]. 

The solid squares, triangles and diamonds in Fig. 0] 
show the total energy E V4 (N) calculated by the DMC 
method for N = 4 trapped particles interacting through 
the i>4 potential with r c — 0.007a^ o as a function of the 
zero-energy scattering length a(0) for a varying number 
of two-body s-wave bound states. The solid squares show 
energies for two-body potentials that support no s-wave 
bound state. For the two-body potentials considered, 
the corresponding four-body system supports no state 
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FIG. 4: (Color online) Solid squares, triangles and diamonds 
show the total energy E V4 (N) calculated by the DMC method 
for four trapped particles interacting through the i>4 potential 
with r c = 0.007ah o as a function of a(Q)/a,ha- The dotted, 
dashed, and dash-dotted lines show the perturbative energies, 
Eq. (|17[) . for two-body potentials that support 0, 1 and 2 
bound states. Open circles show the total energy E a ( )(N) 
calculated by the DMC method for four particles interacting 
through the zero-range pseudo-potential v ps with g = a(0). 
For all cases, the error bars of the energies are smaller than 
the symbol size. The inset shows an enlargement of the region 
around a(0) = 0.03ah o . 



with negative energy, and the hardcore boundary condi- 
tion implied by is met by setting the parameter b in 
Eq. (fT5")l equal to r c . The energies for those V4 potentials 
that support one and two s-wave bound states are shown 
by solid triangles and diamonds, respectively. In these 
cases, the four-particle system supports negative energy 
states, and the parameter b is chosen to coincide with 
the r value at which the free-space zero-energy two-body 
scattering solution has its first and second node, respec- 
tively. This construction of the many-body nodal surface 
assumes that at most two particles scatter at any given 
time, and that the nodal line of the two-body scatter- 
ing solution is not modified by the presence of the other 
atoms [U, [25[ . This "binary approximation" is expected 
to be quite accurate in the low-density regime consid- 
ered throughout this paper. Figure |4] shows that, for a 
given zero-energy scattering length a(0), the energy E Vi 
of the lowest-lying gas-like state increases as the number 
of two-body s-wave bound states, or equivalently 13^/aho, 
increases, similar to the behavior found in Sec. QI] for the 
two particle case [see Fig.rjfc)]. 

For comparison, we consider the energy of N parti- 
cles interacting through an energy-dependent zero-range 
pseudopotential, within first order perturbation theory, 



hu 2 



N(N-1) [2 a(k) 



7T dho 



(17) 



For weakly interacting systems, i.e., for small scatter- 
ing lengths, a{k) can be approximated by Eq. ^ with 
k corresponding to the trap energy scale of 3/2fku. The 
perturbative results for N — 4 particles are shown in 
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Fig. |4] by dotted, dashed and dash-dotted lines for the 
cases when the two-body potential supports zero, one 
and two bound states. The agreement between the per- 
turbative and DMC energies is reasonably good over the 
range of scattering lengths considered. In particular, the 
perturbative expression with the energy-dependent scat- 
tering length predicts the up-shift of the energies with 
increasing number of two-body bound states for a fixed 
a(0) correctly but does not fully capture the change of 
slope of the energy with increasing |a(0)|. 

Next we consider the DMC results for TV trapped 
atoms interacting through the pseudo-potential v ps with 
g = a(0). The potential v ps possesses one bound state 
for g > and no bound state for g < 0. For positive 
scattering lengths, we use the nodal surface of the free- 
space two-body scattering solution for v ps to determine 
the parameters of the Jastrow factor F, Eq. (fl"6)) . i.e., we 
use c = a(0) and d = S(k) with a very small k value. 
In addition to positive a(0), we consider negative a(0). 
Using %\jy with c, d and k chosen so that the boundary 
condition implied by the zero-range pseudo-potential is 
satisfied whenever one of the interparticle distances 
is zero, numerical instabilities associated with large neg- 
ative DMC energies arise. These instabilities are most 
likely associated with the Thomas collapse (2(| H3] , which 
is known to occur for systems with three or more particles 
interacting through v ps with g < 0. 

The open circles in Fig. [4] show the total DMC energy 
_E a ( ) (4) as a function of the zero-energy scattering length 
a(0) for a(0) > 0. We find that £ , a ( )(4) is smaller than or 
equal to E Vi {A) for all a(0). To quantify to which extent 
the Hamiltonian with the shape-independent potential 
reproduces the properties of the Hamiltonian with the 
shape-dependent potential, we consider the energy differ- 
ence AE(N), AE(N) = E Vi (N) - E a{Q) (N). For N = 4 
and 10, we find that AE scales — as might be expected 
for a weakly-interacting Bose gas — with the number of 
pairs, i.e., AE(N)/fou> « N pair Ae V4 , where Ae V4 denotes 
the energy difference introduced in Sec. Ill and N pa i r 
the number of pairs, N pair = N(N - l)/2. For N = 4 
and a(0) = 0.03a/j O , e.g., we find E V4 — 6.1457(2)fiu 
(triangles in the inset of Fig. 4) for the V4 potential 
that supports one two-body s-wave bound state and 
-^a(o) — 6.14189(3)/iu; (open circles in the inset of Fig. 4) 
for the zero-range potential with g = a(0), and thus 
AE = 0.0038(2)^. For comparison, the corresponding 
quantity N pa i r Ae equals O.OOAOhuj. In addition to the V4 
potential, we consider the vq potential. In this case, the 
energy difference AE for comparatively small a(0) is of 
the same order or larger than the statistical uncertainties 
of our DMC energies and, although expected to be valid, 
we cannot explicitly confirm the scaling of AE with N pa i r 
for the vq potential. 



To further understand the implications of the energy- 
dependent scattering length a(E) for N > 2, we deter- 
mine the a(E) that, if used to parametrize the interaction 
strength of the zero-range potential in the many-body 
Hamiltonian, reproduces the energy E Vi . For example, 
to reproduce the four-particle energy E Vi — 6.1457(2) 
for the U4 potential with a(0) = 0.03a^, o that supports 
one bound state, the strength of the pseudopotential 
has to be a(E) = 0.03082(4). For the v A potential, 
this a(E) corresponds to a two-body scattering energy of 
1.53(16)fiu. Thus, the relevant scattering energy for two- 
body collisions occuring in the weakly-interacting many- 
body is, not unexpectedly, approximately given by the 
trap energy scale of 3/2fkv. 



V. CONCLUSION 

This paper studies trapped bosons interacting through 
attractive power-law potentials with r~ 4 and r~ 6 tails. 
For two particles, the energy-dependent pseudo-potential 
accurately reproduces the energies for both shape- 
dependent potentials. Further, we find that the de- 
viations between the energies obtained for the energy- 
independent pseudo-potential and the shape-dependent 
potential scale for small a(0) as (fle/aho) 3 and (P^/aho) 2 
for the potentials with r~ 6 and r~ 4 tails, respectively. Fi- 
nally, we use Monte Carlo methods to extend the treat- 
ment to more than two trapped particles. Again, we 
find that the energy for the shape-dependent power-law 
potential can be reproduced accurately by the energy- 
dependent pseudo-potential if the energy-scale entering 
the pseudo-potential is chosen properly. In general, this 
leads to a self-consistent many-body framework that con- 
siders only binary interactions but includes many-body 
correlations. 

Even if the r -4 results are not directly applicable 
to present-day experiments (combined atom-ion systems 
have not yet been trapped), our comparative study of 
the energetics for r -4 and r~ 6 potentials provides in- 
sights into weakly interacting systems in general and 
van der Waals r~ e interactions in particular. Our cal- 
culations suggest that three-body terms [28| are very 
small in the dilute limit considered throughout this work. 
It seems feasible that the description of systems with 
longer-ranged interactions at the mean-field level can 
be improved by including the energy dependence of the 
scattering length, similar to the frameworks outlined in 
Refs. [lEI- 
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